Nonergodic thermodynamics of disordered ferromagnets and ferroelectrics 
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I. INTRODUCTION 

Many types of disorder in ferromagnets and ferro- 
electrics can induce spin- or dipole-glass states if emer- 
gent random interactions cause frustration and make the 
ground state of a crystal highly degenerate^. Besides 
the classic example of the solid solution (alloy) of fer- 
romagnet with antiferromagneti^ it can be the random 
anisotropySi, simple dilution in the case of competing 
nearest and next-nearest interactions^ or in the presence 
of strong dipole-dipole ones^ as well as structural defects 
such as dislocations- . Now it is widely recognized that 
the unusual nonergodic properties of disordered ferro- 
magnets and ferroelectrics in their glassy phases result 
from the appearance of a number of metastable states 
each having more or less random directions of local spon- 
taneous electric or magnetic dipole moments^. Using the 
different protocols of changing temperature and external 
field (e. g. field-cooled, FC, zero-field cooled, ZFC, etc.) 
one may arrive at the metastable states with different 
net magnetization (polarization) and differences in other 
thermodynamic functions. 

First adequate theoretical description of such history 

- dependence of system properties in glassy phases was 
achieved in Refs. 0, Q using numerical simulations in 
the local mean- field theory of random-bond magnets. It 
was explicitly shown that various regimes of temperature 
and field variations do bring system to different minima 
of thermodynamic potential thus providing the history 

- dependence of thermodynamic functions. Under the 
assumption of the possibility of sufficiently slow (quasi- 
static) variations of external parameters making the pre- 
dictions of static thermodynamic theory valid the number 
of nonergodic features of glassy phases were described in 
Refs. the temperature dependencies of FC and 
ZFC functions, thermal and isothermal remanent mag- 
netizations, the form of hysteresis loops, FC heat capac- 
ity in different fields. These results are in reasonable 
agreement with the experimental data for a series of ran- 
dom magnets. Later the method of the local mean-field 
simulations were used to reveal the presence of noner- 
godic reentrant spin-glass phase in 2d XY random-bond 
model^. 

Subsequent studies of nonergodic phenomena in spin 
systems with random interactions were mainly focused 
on numerical simulations of hysteresis loops. The smaller 



loops inside the main one and evolution via macroscopic 
spin avalanches on the main loop are shown to exist in 
Edwards- Anderson~ and Sherrington- Kirkpatricl i^^d^ 
spin-glass models. The existence of hysteresis loops is 
established in random - anisotropy modet^. The new 
hard-spin mean-field method is developed for frustrated 
random - bond system s -"^"^d^d^ which provides the evi- 
dences for the existence of multiple metastable states and 
can describe a number of nonergodic phenomena with less 
simulation efforts. 

In spite of definite successes in numerical simulations of 
nonergodic effects in glassy phases these methods are not 
destined to provide some general picture for systemati- 
zation and qualitative explanation of these effects. It is 
rather difhcult to find from the existing numerical results 
the possible interrelations between a variety of history - 
dependent phenomena observed in different experimen- 
tal regimes, see Refs. [T7ll25j . Yet it seems that they 
should be necessary present as, e. g., the hysteresis loops 
and ZFC - FC magnetization differences have the com- 
mon origin. Then one may hope to obtain the unified 
picture of such phenomena using analytical calculations 
in some simplified models. But nowadays there are no 
simple enough microscopic models in which nonergodic 
effects in finite field could be described by the analytical 
methods. 

Thus in present situation some phenomenological ap- 
proach may help to achieve general understanding of the 
quasi-static large field response in glassy phases and the 
role played in it by metastable states. Apparently, it 
should be based on some realistic mechanism underlying 
the appearance of such states. Such mechanism for the 
random-interaction systems has been described in Ref. 
f26]. It consists in the subsequent condensations of nu- 
merous sparse fractal modes defined by the eigenvectors 
of the random matrix of pair-wise interactions (exchange 
matrix in magnets or the matrix of harmonic interac- 
tion of polar atomic displacements in ferroelectrics) at 
their localization threshold^i. In random ferroelectrics 
(relaxors) these modes play the role of the notorious 
" soft mode" as they are the delocalized eigenmodes of po- 
lar atomic oscillation with the lowest frequencies. Their 
amplitudes throughout the disordered crystal can be ob- 
tained by the numerical methods of Refs. [l^]- The con- 
densation (freezing) of these modes in relaxors results in 
appearance of stable polar atomic displacements same 
as in the case of ordinary soft mode albeit these dis- 
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placements are highly inhomogeneous and sparse. Their 
analogs in random Ising magnets can be related to the de- 
localized relaxational magnetic modes which freeze first 
at a macroscopic spin glass transition. 

It is important to realize that only the freezing of 
these macroscopic modes can result in emergence of sta- 
ble spontaneous local moments at a glass transition. Such 
objects as nano-domains will never have spontaneous mo- 
ments stable for a macroscopic time so this notion is use- 
less for the description of the true thermodynamic tran- 
sition into the glass phase. Thus Burns temperature in 
relaxors at which experiments reveal the appearance of 
(fluctuating) polar nano-domains marks the onset of the 
paraelectric Griffiths' phase in which only non-analyticity 
in the field dependence of polarization^* and slow non- 
exponential relaxatior^ appear. 

The suggested mechanism of a macroscopic spin glass 
transition is in sharp contrast with that of phase transi- 
tion in homogeneous crystals. The last takes place due to 
the condensation of just one periodic or uniform in space 
eigenmode which corresponds to the largest eigenvalue of 
exchange matrix (in magnets) or the lowest eigenvalue of 
the matrix of effective harmonic interaction (soft mode 
in ferroelectrics). In ideal crystal it is sufficient to sta- 
bilize the other modes with eigenvalues close to the crit- 
ical one throughout some low-symmetry phase regions. 
This stabilization of the rest of near-critical modes re- 
sults from nonlinear mode couplings which are always 
present if there is some finite spatial overlap between the 
critical mode and near-critical ones. As all eigenmodes 
of an ideal crystal is just the plane waves strongly over- 
lapping throughout the whole sample, the condensation 
of just one critical mode makes the system stable. 

In a strongly disordered crystal with random interac- 
tions the situation changes principally. The randomness 
of pair-wise interactions makes the eigenmodes near the 
boundary of the spectrum localized^! so the transition 
may start only with the condensation of the first de- 
localized eigenmode at the localization threshold some- 
what far from the lower boundary of the spectrum. Ow- 
ing to the sparse fractal structure of the modes near 
this threshold^l, the condensation of the first delocal- 
ized mode can not stabilize the other sparse near-critical 
modes which do not overlap essentially with the first one. 
So they proceed to condense at slightly lower tempera- 
tures until almost all crystal sites acquire spontaneous 
moments. If the average fractal dimension of these con- 
densing modes is d/ < d , the average number of sites 
participating in each mode is of the order iVi = N'^f^'^. 
So the number of them needed to cover the iV-site crystal 
is iVo = N/Ni = N^-(df/d)^ 

Thus strong disorder can transform the ordinary phase 
transition into the macroscopic number of them spread- 
ing in some temperature interval. As each mode after 
condensation can have (at least) two stable states with 
the reverse directions of local moments, we may end up 
with up to 2^" metastable states. Note that each of 
these numerous subsequent transitions is still a macro- 



scopic one as the macroscopic number of local sponta- 
neous moments of order A^i appears at it. Apparently, 
the mean-field estimate of the potential barrier between 
metastable states would also give the macroscopic value 
proportional to iVi. 

Here we may note that above picture describes only the 
possibility of emergence of numerous metastable states in 
crystal with random interactions. Do it actually realizes 
depends on the details of pair-wise interaction between 
the sparse fractal modes. If this interaction tends to ori- 
ent the net moments of the modes in parallel, we may 
have ordinary ferromagnetic (ferroelectric) phase with or 
without some metastable states. Principally, sketched 
above mechanism allows for phenomenological descrip- 
tion of the thermodynamics of metastable states in the 
spirit of Landau's theory of phase transitions. The case 
with purely glass transition was considered in Ref. [2^ . 
Here we present the simplified and more general deriva- 
tion of the phenomenological Landau's potential for the 
generic case of competing ferromagnetic (ferroelectric) 
and glassy interactions in disordered crystals. For a large 
class of such random systems it is possible to obtain semi- 
analytical description of the properties of the emergent 
metastable states thus providing the comprehensive ac- 
count of possible nonergodic effects in random ferromag- 
nets and ferroelectric relaxors. 



II. LANDAU'S POTENTIAL FOR RANDOM 
FERROMAGNETS AND FERROELECTRICS 

Further we mainly resort to the magnetic terminol- 
ogy and designations just to be specific. Let us consider 
the ferromagnet with the second-order transition and 
one-component order parameter where some disorder of 
random-interaction type presents tending to destroy the 
ferromagnetic order. According to above considerations 
we can introduce in this case the Landau potential for 
some specific realization of disorder depending on the net 
magnetizations of sparse fractal modes, TOj, i ~ 1, A^o- 
If the temperature interval (Tg , Tg — AT) at which 
the condensation of these modes takes place is narrow, 
AT <C Tg then we can consider the small m.i near Tg and 
expand the potential F (m) in powers of m^. It would 
contain only even powers of these order parameters due 
to the global inversion symmetry. Also we can retain in it 
the terms with no more than two different modes as one 
can suggest that their sparse spatial structure makes the 
simultaneous interactions of three or more modes negli- 
gibly small. 

Further we may note that AT <C Tg is actually the nec- 
essary condition when one consider the frustrating disor- 
der causing the (near) degeneracy of condensing types 
of magnetic (dipole) arrangements. This means that in 
such case the condensing modes have nearly the same 
thermodynamic potentials, close transition temperatures 
and thermodynamic parameters. The simplest way to 
imitate this near degeneracy is to suppose that all ther- 
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modynamic properties of all modes are identical. Then 
F (m) should be symmetric under all permutations of m^. 
The immediate effect of this assumption will be the merg- 
ing of actually subsequent transitions to just one point, 
Tg. In this way we avoid the consideration of narrow 
temperature interval {Tg , Tg — AT) and gain the sig- 
nificant simplification of the model which still preserve 
the essential features of nonergodic transitions caused by 
frustrating disorder. We may say that model with per- 
mutation symmetry is the "minimal" one and, probably, 
the simplest model allowing for more or less adequate ac- 
count of real nonergodic effects in large class of random 
ferromagnets and ferroelectric relaxors. 

At last, the quadratic terms of the potential should 
refiect the competition between ferromagnetic (ferroelec- 
tric) and glassy orders. Thus we arrive at the unique 
form of potential up to forth order in rrii obeying the 
above criteria, 

/ TVo \ ^ No 

\z— 1 / — 1 

No A / \ ^ 

+iE-^ + i^ E-n -^T.m. (1) 

Here h is external field conjugate to the net magnetiza- 
tion of a sample 



the same composition. Here we should note that this is 
true only for the solid solutions with deeply frozen disor- 
der. In the relaxors with annealing mediated ordering of 
components such as PSN the parameters of the potential 
will also depend on the degree of ordering achieved at the 
high-temperature annealing. 

Generally, the crystal with negligible disorder should 
have Tg <C Tf as in this case we may have just the or- 
dinary transition into homogeneous ferro-phase without 
any traces of disorder. So the growth of disorder will 
result in the increase of Tg along with the lowering of Tf . 

The difference between Tf and Tg must be small to 
make the expansion of F in small rui meaningful. Ac- 
cording to Eq. ([T]) below Tf ferromagnetic order may set 
in, while below Tg glassy states described by the iVg — 1 
- dimensional order parameter composed of the indepen- 
dent TTii — rrij components can appear favoring the an- 
tiparallel m; orientations. What actually results from 
the competition of glass- and ferro-order depends on the 
relation between Tf and Tg values as well as the rates 
of Tf and Tg decreasing. Finding the evolution of F (m) 
minima for different Tf and Tg relations we can get the 
possible variants of phase sequences in crystals with ran- 
dom interactions. In spite of huge amount of glassy min- 
ima that F (m) may have, to find all of them is rather 
simple task in the present model which can be fulfilled 
by semi-analytical methods. 



No 



Note that this relation is the consequence of the adopted 
principle of mode equivalence which includes the supposi- 
tion of equal numbers of sites participating in each mode. 
Nq in the denominators of Eq. (1) ensure the equal order 
of different terms' contributions at large Nq ^ 1. Note 
also the absence in (1) the term 
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-T 



TO, 



No 



allowed formally by the above criteria. It can be removed 
by the linear transform of resulting in rescaling of 
h and Tf and appearance of two forth-order terms with 
three- and four-mode interactions which must be dropped 
due to assumed smallness of such interactions. Thus the 
presence in F of this term amounts just to the h and Tf 
rescaling and it can be omitted. 

The coefficients a and b in F are some constants spe- 
cific for a given disorder realization, while Tf and Tg are 
linear decreasing functions of temperature T changing 
their signs at temperatures Tf and Tg correspondingly 
also being disorder dependent. Yet in the ferroelectric 
and ferromagnetic solid solutions all potential parame- 
ters must be the self- averaging quantities. It means that 
they depend only on the impurities' concentrations in ac- 
cordance with the experimental data showing no notice- 
able variations of the properties of different samples with 



III. THERMODYNAMICS OF COMPETING 
GLASS- AND FERRO-STATES IN ZERO FIELD 



_ No 

Using the notations [to*'] ~ Nq ^ to*', we can rep- 

1=1 

resent the potential density /(m) = F{m)/No in the 
simple form 



9 r^21 , Tf Tg^^ ^ a ^^4^ ^ b ^^2^2 _ ^ _^ 



-TO , 
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Differentiating / (m) with respect to rrii we get the equa- 
tions of state 

[Tg + b [to^] ) TOi + am^ = h + {Tg — Tf) rh (2) 

The solutions of these equations are minima of / (m) 
describing possible (meta)stable states if they render the 
positive definiteness of the matrix 



^ gV(m) 
dmidrrij 



-f A^o"^ {Tf — Tg + 26TOiTOj) . 



Let us consider first the homogeneously magnetized 
ferro-state with the equal to; = too. From Eq. ([2]) we 
have 



T/TOo + (a + 6) TOq = h 



(3) 
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Thus ai h = the ferro-state appears at r/ < if a + 6 > 
0. It is stable for 



3a + 6 



(4) 



Introducing the variable 



Turning to the possible glassy solutions of Eq. ^ with 
unequal m^, we note that all rrii obey the same equation 
so such solutions exist if there are several real roots to Eq. 
(121). Being of third order with respect to rrii it can have 
one or three real roots. In the last case rrii can acquire 
only two of the root values with the largest modules: one 
positive, m+ > 0, and one negative, m_ < 0, as only 
these two can make the coefficient at 5i_j in matrix G 
positive which is necessary for the stability. Thus every 
glassy state is defined by the number, N^, of rrii having 
values (or the number of m_ ones, — Nq — N^). 
Hence in the glassy states 



n± = N±/No, 



m - 
= 1 



n_|_r7i_|_ + n^rri- 



and to find a glass state with a given (ri_) we need 
to obtain the corresponding roots, m±, of Eq. (2) which 
now becomes 



(rg + b [m^] ) m± 



(5) 



Thus the present model may generally have the quasi- 
continuous set of states (local minima) defined by the 
parameter n+ which changes in the interval < n_|_ < 1 
in the infinitesimal steps ±Nq^. Note that there are 
No 



,j , states with the same iV+ n+-/Vo which differ 

by the permutations of rrii and have the identical ther- 
modynamic parameters due to the adopted permutation 
symmetry. 
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we get from Eq. (5) the following equations 



1 ] (n+ — xri-) R (x, n+, /?) 

^9 



X {x — 1) i? (x, n+, H — ( 

a \ -To 



3/2 



(6) 



m+ = ^-Tg/aR {x,n+, /3) (7) 
R {x, 71+, (3)= [l-x + x^ + (3 (n+ + "^^^ 

Here /3 = b/a. 

Thus we need to solve only one Eq. (6) for x — 
X [n-\-,Tg,Tf ,h) to obtain using Eq. (7) and x defini- 
tion full description of the thermodynamic properties of 
a crystal in the metastable states with a given . So we 
have for the net spontaneous moment in such states 



m = (n_|_ — xri-) ^—TgjaR (x, n+, /?) (8) 
and for the Edwards- Anderson glass order parameter 

q = [to^] — TO^ = {~Tg/a) n^n_ (1 + x)^ R {x, n+, . 



Further we consider the simplest case a > 0, 6 > so the 
glassy metastable states appear at < . 

Then we have for the susceptibility X = ^ 



X ^ ^ Tf - Tg - TgR{x,n+,/3y 



(1 -I- x) {2x - 1) (2 - x) + 2/3 [(2a; - 1) n+ + x^ (2 - x) n_] 
{2x - 1) n+ + (2 - x) n_ + 2/3n+n_ (1 + x) 



(9) 
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The equilibrium value of the potential for a given state 



feq {n+, Tg,Tf, h) = [Tgq -f T/TO^ - ihwi) /4 

Here we omit the term arising from the degeneracy of 
glassy states 



T 

'n' 



"^config 



N \ N+ 



No 
N 



as it vanishes at large N. Then we get for a given 
metastable state the entropy 



o _ df^g{n+.Tg,Tf,h) 



\{T'gq + r'j, 



'l,g- dT 



and the heat capacity 
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C , I i\i -} I / I i\ mm+ \n+ (2x — 1) — ri-x (2 



TgX ^' n+{2x-l) + n^{2-x) + 2n+n^f3{l + x) 

ml [n+ {2x - 1) + n_a:^ (2 - x)] + n+n. (1 + x) g-^ (t/ ~ Tg) 
^^3' n+{2x-l) + n^{2-x) + 2n+n-f3(l + x) ^ ' 



For a > 0, 6 > the stability conditions for the so- 
lutions of Eq. ^ providing the positive definetness of 
matrix G are represented by the inequalities 



1/2 < a; < 2, x > 



(11) 



Note that the condition x > ensures the stability with 
respect to the ferromagnetic fluctuations while at the 
boundaries of the first inequality the glassy instability 
occurs. 

Turning to the case /i = we can see that simple solu- 
tions of the equation of state (6) exist for n+ ^1,0 and 
n+ = n_ = 1/2. In the first case Eq. ([6]) tends to the 
zero-field limit of Eq. Q and we have 



5^-r;.mV2,C=^ 



and Eqs. lfTTj) defines the sector on the (ry, rg)-planc in 
which nearly homogeneous states with n-^- —>■ 1,0 are 
stable 



3 + P 3 + 4/3 



(12) 



Note that such states exist only in the part of the region 
Eq. ([4]) where the homogeneously magnetized ferro-state 
is stable. This just means that the set of metastable 
states to which the states with n+ 1,0 belong has 
more narrow region of existence. 

For the state with fully disordered moments = 
n_ = 1/2 ) we have x = 1 from Eq. (I6|) at ft, = so 

^ = -'^g9/2, C = 2{a+b) 



and X > requires that 



3 + /3 



(13) 



along with < 0. It can be shown that the sector of 
stability of fully disordered states is the largest one and 
the stability regions of all other states belong to it, the 
most narrow sector being that of nearly homogeneous 
states ([T2|) . 

Further we find 



9/, 



eq 



dnj 



(m_|_ + TO_) (to+ — m_)'^ . 



Thus /eg (n+ ,Tg,Tf,h) as function of 
tremum at 



has one ex- 



m+ + m_ = 



and we have at this point 



(x=l), 



eq 



8aTO+x(T-/ - Tg) 



1 + /3 



1 + Afin+n^ 



So it is minimum at Tg < Tf and maximum at Tg > Tf. 
Thus the states with ri+ = n_ = 1/2 have the lowest 
potential at Tg < Tf while the global minimum is at the 
states with n+ —^ 1,0 at Tg > Tf. 

Hence the first-order transition from the fully disor- 
dered glass to ferro-state takes place at Tg — Tf. One 
can see from the above equations that potentials of these 
states indeed become equal at Tg — Tf. The phase dia- 
gram of the model is shown in Fig. [T] The topology of 
the phase stability regions in it is the same that is found 
in the various microscopic models of disorderd magnets, 
cf. Refs. [H, [iOlsl. 

In a specific crystal there is some linear relation be- 
tween Tf and Tg both being linear functions of T. It 
defines a straight line on (Tf,Tg) plane which crystal fol- 
lows under T variations and, hence, the phase sequence 
proper for a given crystal. For example, if the relation 
Tf = Tg + To holds, To being some constant, the crystal 
undergoes just one second-order transition - into the ordi- 
nary ferro-phase if tq < or in the glass phase if tq > 0. 
In the first case there is another qualitative change of 
thermodynamic properties (not a phase transition) - the 
subsequent appearance of numerous metastable states in 
addition to ferro-state begins when crystal reaches the re- 
gions of their stability belonging to the sector p3|) . Yet 
in these cases the application of field causes no drastic 
changes in the temperature dependencies of ffi, x and 
C apart from the usual differences in their EC and ZFC 
values. 

Much more pronounced field effects occur when the 
path on (Tg,Tf) plane intersects the line Tg = Tf. Then 
glass or ferro-phase can appear as intermediate one in 
some finite temperature interval neighboring para-phase 
and first-order transition between them takes place at 
Tg — Tf. Yet there are also different temperature behav- 
ior when path crosses the region of coexistence of ferro- 
and glass states (the sector in Fig. 1 bounded by the 
dashed lines) and when it always stays in this region. 
These four principally different paths are shown in Fig. 
1. Dotted lines in it correspond to paths Tf = 2.5x^-1- 1.5a 



/ para-phase 




/ para-phase 




FIG. 1: (Color online) Phase diagram of the model for /3 = 5 (a) and /3 — 0.1 (b). Thick lines denote the phase transitions 
between the phases. Upper dashed lines show the stability boundaries for the ferro-state (Eq. |4| and lower dashed lines show 
that of fully disordered states with n+ = 0.5 (Eq. I13|l . Dotted lines correspond to paths Tf = 2.5rg + 1.5a (A and Ai) and 
Tf = 0.5(rg - a) (B and Bi). 



(A and Ai) and r/ — 0.5(Tg — a) (B and Bi) which we 
use further to present some typical temperature varia- 
tions of thermodynamic variables which differ essentially 
on them. 

The spin-glass phase on paths B and Bi is often called 
"reentrant" as crystal returns again on cooling to the 
equilibrium state with zero net magnetization. The pos- 
sible microscopic mechanism of such reentrance is con- 
sidered in Ref. 0. 

We must note that here the first-order transition does 
not mean that at Tg = Tf the crystal do jump from one 
state to another on the laboratory time scale. Actually 
it must traverse the rough potential landscape and to 
overcome on the way a number of macroscopic barriers 
between the local minima (partially ordered states) to ar- 
rive to the global minimum. It takes an infinite time and 
the standard thermodynamics indeed predicts formally 
what will happen after infinite time when system will be 
fully relaxed. But actually one has just finite laboratory 
time to measure the thermodynamic quantities. 

It is quite usual situations in all nonergodic systems (i. 
e. those with metastable states) - they all have infinite 
(Arrhenius) relaxation time for the decay of metastable 
states into stable one irrespective of their number - to tra- 
verse just one macroscopic barrier also needs an infinite 
time the same as to traverse many of them. This cir- 
cumstance causes no difficulties in the application of the 
results of the equilibrium thermodynamics to the known 
nonergodic systems such as ordinary ferroelectric and fer- 
romagnets in their ferro-phases. The simple example is 
the description of hysteresis loop. Here the infinite-time 
thermodynamic prediction that spontaneous magnetiza- 
tion should reverse its direction just when the field change 
its sign is easily modified with due account of the long- 
living metastable states (in which field and magnetization 
have opposite signs) . As they can not relax to the global 
minimum at the laboratory time and persist until they 
become unstable (at the coercive fields) one just include 



them in the description of the response to sufficiently 
slow (quasi-static) field thus obtaining upper and lower 
branches of hysteresis loop. Apparently, to apply the re- 
sults of static thermodynamic to the dynamic process of 
field variation we must be sure that during it the sys- 
tem stays close to a local minimum. Hence, in this case 
the quasi-static condition means that the characteristic 
time of field variations is greater than the magnetization 
relaxation time in a local minimum. 

In the same way we can get the description of quasi- 
static responses from the thermodynamic results for the 
present model. To do this one must assume that quasi- 
static regimes can be attained experimentally, i. e. that 
sufficiently slow temperature or field variations can be 
achieved on the laboratory time scale which ensure the 
location of system near some minimum during these vari- 
ations. Then we just should take into account that once 
the system entered the metastable state it can not leave 
it in such quasi-static processes until that state becomes 
unstable. 

Hence, the validity of this picture of the temperature 
evolution in which crystal is trapped in a local minimum 
and has the obtained above thermodynamic parameters 
crucially depends on the rate of cooling or heating. It 
should be sufficiently small to provide the quasi-static 
evolution of the system, that is the characteristic time of 

temperature variations, to = ( '^^"^ ) : must be greater 
than the largest magnetic relaxation time in the specific 
minimum. The same is true for characteristic time of 
field variations. One may expect that when we are not 
close to Tg nearly homogeneously magnetized states with 
TO ~ Too would have rather moderate relaxation times 
comparable to those the crystal has deep in para-phase. 
Also far from Tg the magnetic susceptibility of highly 
disordered states with to <C toq may be very small mak- 
ing the amplitude of magnetization relaxation in them 
small too. Then we would observe their nearly static pa- 
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rameters even for to less than relaxation times which can 
spread up to huge macroscopic values in such states. If 
X reaches large values close to Tg in the states with low 
rh <C ruQ their thermodynamic description would become 
useless due to the strong aging effects in them. Then one 
should necessarily resort to the dynamic theory to get 
the adequate account of their (dynamic) properties. 

Provided the above quasi-static conditons are fuUfiUed 
another question inherent to such highly nonergodic sys- 
tems arises - to which of numerous stable states the sys- 
tem will go after its present state becomes unstable? 
Strictly speaking, one should again turn to the dynamics 
to answer it, yet, as we will see, there is actually no vast 
choice of possibilities in the framework of the heuristic 
assumption that system adopts the most smooth varia- 
tions of its average magnetization. 

Thus crystal can stay in metastable state for a very 
large (Arrhenius) time and leave it only when it becomes 
unstable. Before this the only way to bring the noner- 
godic crystal to some other state is to apply the external 
field which could make the present state unstable. This 
circumstance made former investigators to believe that 
ferroelectric relaxors have no phase transitions exhibit- 
ing only "field - induced" ferroelectricityil. 

All these considerations must be taken into account 
in the interpretation of the temperature dependences of 
metastable parameters described by the above equations 
and shown in Figs. - O Note that to plot them one 
does not need to solve the equation of state ^ as it 
provides along with the expression for a thermodynamic 
variable the parametric representation of its temperature 
dependence, x being the parameter varying in the region 
of stability of a given state. 

The temperature dependencies of metastable states' 
magnetizations in Figs. [5] - [5] are bounded by the lines 
at which states become unstable. In the glass phase they 
correspond to x = 1/2,2 so their equations follow from 
Eqs. [ni[5]at these values of x 

Here minus and plus signs are for (Figs. [21 H]) and (Fig. 
121) correspondingly. In the ferro-phase this relation is not 
generally valid as some states seize to exist due to the 
ferromagnetic instability (x — > oo). In this case it is 
hard to get the compact analog of Eq. (fT4l) 

As crystal leaves its state only when it becomes unsta- 
ble, it is rather evident that the quasi-static heating (Fig. 
121) or coohng (Fig. ^ of crystal being in the ferro-state 
beyond its stability point in the glass phase makes it to 
go through the succession of metastable states on the 
boundary of their stability regions. So the self-organized 
criticality shows up in these processes as system is per- 
manently unstable in finite temperature intervals and its 



evolution proceeds via small but macroscopic avalanches 
of spin upturns. 

Eqs. (fT4l) describe the temperature dependencies of 
magnetization in these regimes. In particular, for the 
vanishing of magnetization on heating in Figs. |21 HI we 
get at Tg ^ 

\Tf\\3a + 2bJ ■ 

Also the specific temperature hysteresis occurs in the 
cooling-heating-cooling cycles (Fig. [21) and cooling- 
heating ones (Fig. ^ starting on the paths A and B from 
the high-temperature phases (in the presence of infinitesi- 
mal positive field to have m > 0) . These cycles are shown 
by the directed lines on the temperature dependencies 
of thermodynamic parameters. Their remarkable prop- 
erty is that they allow to enter a variety of metastable 
states (actually all in the range 0.5 < n-|- < 1 in the 
case of intermediate glass phase in Fig. [2] and those with 
0.628 < n+ < 1 having the limited stability regions in 
Fig. [21). For example, in the case of Fig. [2] we can stop 
heating at some T when crystal traverses the stablility 
boundaries of metastable states and start cooling it thus 
trapping crystal in the state it was in at the stop. 

Note the sharp rise of the susceptibilities and heat ca- 
pacities of metastable states in Figs. [21 [SI Actually this 
denotes the divergence of x and C (cf. Eq. fTOl) ) due to the 
ferromagnetic instability metastable states experience in 
the ferro-phase at Tf > Tg {Tg/a < —1 in Fig. [21 and 
Tg/a > -1 in Fig. [21). 

On the path Ai the appearance of intermediate glass 
phase in zero field is marked by the jumps in heat capac- 
ity and X temperature derivative (Fig. [JJ while on the 
path Bi the intermediate ferro-phase (Fig. has usual 
Curie - Weiss anomaly of x and the usual jump in heat 
capacity (not shown). In these cases the appearance of 
metastable states can not be observed in zero-field experi- 
ments - at all T the crystal will stay in the state it enters 
first on coohng. But the presence of metastable states 
can be easily revealed in the glass phase. Here stopping 
at some T and applying for some time the field (above 
some threshold one needed to make the former state un- 
stable) one can end up in every metastable state choosing 
the value of the applied field. In the further temperature 
variations the crystal will stay in the so prepared state 
until it vanishes. 



IV. THERMODYNAMICS IN FINITE FIELD 

Further we examine the properties of metastable states 
in finite fields. In the glass phase at Tg < Tf all 
metastable states always have x > so they have not 
ferromagnetic instabilities here. Again, substituting x — 
1/2 and a; = 2 in Eqs. ([SI [HI) we get the parametric 
representation of two lines on (m, /i)-plane bounding the 
region where these states exist (0 < n+ < 1 being the 
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FIG. 2: (Color online) Temperature dependencies of metastable states' magnetization, susceptibility and heat capacity in zero 
field for /3 = 5 corresponding to path A in Fig[lla). Dotted lines show the stability boundaries and dashed ones represent the 
parameter's variations predicted by the equilibrium thermodynamics (those of the states with n+ — 1/2 in the present case). 
Directed lines show the evolution of magnetization in real-time quasi-static regime. 




-2 Xg/a 



FIG. 3: (Color online) The same as in Fig. [2] for path B in Fig. \Tla) 




FIG. 4: (Color online) The same as in Fig. [2] for path Ai in Fig[llb). Dash-dotted line corresponds to C/T. x s-nd C are 
shown only for the n+ = 1/2 state as others are inaccessible ai h — 
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Xg/a 



FIG. 5: (Color online) Temperature dependencies of 
metastable states' magnetization in zero field for /3 = 0.1 
corresponding to path Bi in Fig. Hlb). Dashed directed line 
show the decay of nearly homogeneous state's magnetization 
in the infinitesimal field. The other lines have the same mean- 
ing as in Fig|21 



parameter). Excluding n+ we get the equation two real 
roots of which define the couple of these lines 



3a + 6 



2a 



2a 



Tf)m 



1/3 



2/3 



= 



(15) 



These roots have h + {rg — Tf) fh > and h + 
{Tg — Tf)rh < with the interval < rt+ < 1 be- 
ing swept when h varies correspondingly in the intervals 
—hf < h < he ,—he < h < hf 



hf 



3a 6 



3a + 46 



4r. 



1 + 13 
'3 + /? 

1 + /3 
'3-^4^ 



(16) 
(17) 



The magnetization curves m = fh{h^n^) of metastable 
states defined by Eqs. ^ [S]) fill the region bounded by 
the hne segments of Eq. p2|) as Figs. [H [7] show. Again 
the self-organized criticality appears in these segments 
as system is permanently unstable in them and its evo- 
lution proceeds via the avalanches of spin upturns. Such 
avalanches and small magnetization jumps are shown 
to exist in hysteresis loops of Sherrington-Kirkpatrick 
modeliifi^. 

At \h\ > he all metastable states vanish leaving the 
only stable homogeneous state with magnetization curve 
mo{h) given by Eq. As it was discussed above 

the crystal will stay in a local minimum until it van- 
ishes so Fig. [6] represents possible forms of hysteresis 
loops for quasi-static periodic field. If amplitude of such 
field \h\ > he the form of quasi-static loop is completely 
described by Eqs. [T^ . For smaller amplitudes the 



inner loops inside the outer one will be observed formed 
by the boundaries of Eq. (|12p and fh — fri{h, n+) and 
fh — rh{h,l — n-)_) curves with some n+ defined by the 
field amplitude. 

We can find the analog of the coercive field he for the 
glassy loops putting m = in Eq. ill 



2a 



3a + 2b 



3/2 



The inclined hysteresis loop, see Figs. [HUbjC), El^a), 
is the hallmark of glass phase; such loops are ubiqui- 
tous in random magnets and relaxors. They are found 
in the numerical simulation of spin glass phases of 2d 
short-range random-bond Ising modeU^, Sherrington- 
Kirkpatrick modelilfi^, random anisotropy modeU^. as 
well as in the local mean-field simulations^,—. Qualita- 
tively the transition at tj = Tg manifest itself by the 
change of loop form, inclined in the glass phase it be- 
comes the ordinary one with the vertical sides in the 
ferro-phase as disordered metastable states can not be 
reached in the quasi-static periodic field, see Fig. [D^a). 
It is shown that such change of hysteresis loop takes place 
in Sherrington-Kirkpatrick model with the growth of fer- 
romagnetic exchange^"'^. 

Also smaller loops vanish in the ferro-phase as the field 
evolution of all metastable states shown on Fig. [HJa), 
[7t^b,c) ultimately ends by the transition into the homo- 
geneous one. For & = it happens strictly at Tf — Tg 
while for finite 6 > inclined loop transforms into rect- 
angular one already in the glass phase slightly before the 
thermodynamic transition into the ferro-phase. Gener- 
ally this change takes place when he becomes less than 
—hf, that is at 



l + (3 
3 + 13 



' 3 + 13 
(3 + 2l3f 



It is also easy to get for all h the equilibrium magne- 
tization nieq corresponding to the states with the lowest 
potential having a; = 1 in glass phase as was shown above. 
From Eqs. (O, ([5]) we have 



for the equilibrium states with 



(18) 




The Almeida-Thouless field hAT marks the thermody- 
namic transition between the completely ordered equi- 
librium phase at \h\ > Hat to the sequence of disordered 
ones at \h\ < hAT- 
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FIG. 6: (Color online) Field dependencies of metastable states' magnetizations for (3 — 0.1 on the path Ai, jj. = m\J —aj Tg, 
H = h^^—a/T^, (a) - Tg = — l.lo, (b) - Tg = —0.75a, (c) - Tg = —0.5a 




-0.15 H 0.15 -0.2 H 0.2 -0.3 H 



FIG. 7: (Color online) The same as in Fig. 6 for /3 = 5 and path B. (a) - Tg = —1.5a, (b) - Tg — —0.65a, (c) - Tg — —0.5a 



For the disordered states rheq is represented in the in- 
teriors of hysteresis loops in Figs. (O [7|) by dashed lines 
with tangents defined by the equilibrium susceptibility 



Xe 



dh 



1 



(19) 



and ending dX h = zLHat- 

Generally Xeq is unobservable quantity as to follow the 
relation in Eq. (fT8|) under the field variations the crystal 
has to overcome the macroscopic barriers in series of first- 
order phase transitions (stepping in them from the state 
with some to that with 7i-|_ ± Nq^) which are not 
shown in Figs. [H [7] (see Fig.l in Ref. [IB)]- As we 
discussed above, this would need a very long time and 
crystal will stay in the local minimum on the laboratory 
time scales. So the only susceptibility measurable under 
field variations is that given by Eq. ^ which defines the 
tangents of the magnetization curves ffi = fh{h, n+) of 
the specific metastable state at a given field. It is always 
less than Xeqi thus in the equilibrium states with x ~ 1 



X 



1 



h\j, 



while at the right {x — 1/2) and left (x = 
of hysteresis loop in glass phase we have 

6/3 



2) boundaries 



X 



1 



(3 + 4n±/3) [3 + p{l + 3n±)] 



with plus and minus sign correspondingly. 

Yet in the case b a (/3 ^ 1) Xeq can be determined 
just from the shape of hysteresis loop. Indeed, Eq. (fT5|) 
gives then for its boundaries the almost straight lines 
with tangent equal to Xeq, see Figs. ID^b, c). But for 
large /3 this relation is lost owing to the intricate form 
of the loop as Figs. [Tlja) show. For small (3 the shape 
of hysteresis loop in the glass phase is mostly defined 
by the fields he and hf. In the ferro-phase { Tf < Tg) 
the loop boundaries are no longer given by Eq. psp as 
here some states vanish due to ferromagnetic instability 
before reaching the glassy one at a; = 1/2 or a; = 2. Also 
in the ferro-phase the expression for the field hf, Eq. ([T| 
changes to 



3 V 3 (a -I- 5) 



(20) 



in the interval 



3 + f3 3 + 4/3 

< ^9 < ion I m ^/- 



3(l + /3) 



12(l + /3) 



The field in Eq. is just the coercive field for the 

homogeneous ferro-state. Note that at 



3-^4^ 
12(l + /3) 
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FIG. 8: (Color online) Temperature dependencies of charac- 
teristic fields; (a) - /3 = 0.1 on the path Ai, dash-dotted line 
represents h = 0.2a, (b) - P — 5 and path B. Dashed lines 
- Hat, dotted ones - he- Solid lines show the modules of hf 
and he- 



the field hj and he merge at this value so at 

3 + 



^9 > 



12(1 + 13) 



disordered metastable states no longer exist at all h and 
we have the ordinary ferromagnetic loop. The temper- 
ature dependencies of the loop's characteristic fields for 
paths Ai and B are shown in Fig. [51 in the ferro-phase 
he and he were found numerically. 

Now we can turn to the temperature dependencies of 
metastable states' parameteres in a constant field. In 
general these states gradually vanish at larger fields, first 
the most disordered ones, as one may expect. Fig. [9] 
illustrates this process for the path A. Temperature de- 
pendencies of fh, X and C in constant field are closely 
related to that of hf and he- Thus for paths A and 
Ai hf first grows on cooling and then diminishes after 
reaching the maximum, see Fig. [Slja) for path Ai where 

IT- f, max 

« 0.28a is attained at Tg « -0.24a. Also he > hf 
exhibits the similar behavior. This results in the essen- 
tial changes in x and C dependencies for field-cooling 
and field-heating in Fig. [TO] and, especially, in Fig. [TT] 
as compared to those in Figs. [H [3| for h — 0. In the 
insets of Figs. [TOl [U there are slight anomalies of these 
parameters when crystal enters on cooling the interval 
where h < hf and leaves the homogeneous state to join 
the sequence of glassy ones with n+ < 1 until it reaches 
the first state being stable throughout all this interval. 



It varies for different paths and fields. In Figs. [TOl [TT] 
they are those with n+ = 0.916 and n+ = 0.921 corre- 
spondingly. This process is nothing else than the descend 
throughout the upper branch of the growing hysteresis 
loop. In further cooling the crystal stays in this state un- 
til it vanishes due to ferromagnetic instability near the 
point where h ^ he (at another branch of the loop having 
now the rectangular form, see Fig. [D^a)) for path Ai). So 
X and C diverge at this point according to Curie- Weiss 
law and crystal returns to the ferro-state. This is new 
feature for path Ai which is absent in zero field. Fig. [J] 
On heating from this state again starts diminishing at 
the lower boundary of the interval where h < hf reaching 
n+ = 0.916(0.921) (upper branch of the loop). Further 
heating in this state ends up at the upper boundary of 
h < hf interval in the sequence of states with n+ growing 
up to 1 (lower branch of the loop). This cause the slight 
high-temperature spikes in x and C manifesting the glass 
instabilities of these states. 

Here we may note that Fig. [TT] explains the origin 
of additional pronounced peak in x appearing in relaxor 
PMN in FC regim e^^i^^ . It is the consequence of the 
instability of nearly homogeneous state the crystal en- 
ters on field-cooling. Moreover, the behavior of PMN 
dielectric susceptibility observed in all field-cooling-field- 
heating cyclei^ii^ resembles qualitatively that of x in 
Fig. [TT] The appearance of additional spikes of 
magnetic susceptibility in FC regime was also regis- 
tered in random ferromagnets PrNiQ ^Coo.rO^^ and 

{Feo.l7Nio.83)75PwB6Ak^\ 

We should also note that magnetizations in the regimes 
of cooling from the para-phase shown in Figs. llOf a). 
[TTT a) are the field-cooled ones, mpc- These figures rep- 
resent also the zero-field-cooled magnetizations, mzFC, 
obtained at turning on the field after cooling the sam- 
ple in zero field. They are the metastable curves with 
n+ = 1/2 (if they exist) or the stability boundaries cor- 
responding to the lower branch of inclined hysteresis loop 
or to the homogeneous states lines when the loop is rect- 
angular as in Figs. [TOT a). Indeed, in zero- field cooling 
process on paths A and Ai the system is trapped in the 
ri+ = 1/2 state, see Figs. [Ifa), [l][a), and stays in it af- 
ter turning on the field or goes to the lower loop branch 
or to the ~ 1 state depending on the loop form and 
field strength. So in Fig. [TOTa) mzFC corresponds to the 



lower stability line while in Fig. [TTTa) it follows i 



1/2 



curve and then jump to n+ — 1 curve when former seized 
to exist. 

Here we should note that often in experiments and sim- 
ulations another definition of mzpc is used, namely, the 
magnetization registered in field-heating after zero-field- 
cooling process. As Fig. [TOTa) shows this is the ill-defined 
quantity as it will depend in this case on the temperature 
at which the turn from cooling to heating takes place. Yet 
when the field at the turning point leaves the n+ = 1/2 
state stable as for path Ai in Fig. [TTT a) the field-heating 
after zero- field-cooling protocol gives a unique mzFC the 
same as in simple ZFC process described above. The de- 
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FIG. 10: (Color online) Temperature dependencies of metastable states' magnetizations, susceptibility and heat capacity in a 
field h — 0.13a for P — 5 corresponding to the path A. Insets show the high-temperature anomalies. 



tailed picture of mzFC and mpc behavior for path Ai 
is shown in Fig ll2l The qualitativly similar temperature 
dependences of these magnetizations are found in local 
mean-field simulations^ and in field-heating after zero- 
field-cooling experiments in PMN^®. 

The applied field cause also some peculiarities in tem- 
perature dependencies of jfi, x a-^d C when ferro-phase 
is intermediate between para- and glassy ones (path B). 
As Fig. [13] shows they appear in field-heating regime of 
crystal being initially in the most disordered state with 
n+ = 1/2. After reaching the limit of stability of this 
state the crystal enters the sequence of states with larger 
n+ on the boundaries of their existence regions until it 
joins the fully homogeneous ferro-state with n_|_ ~ 1. 
This process is accompanied by the spikes in x and C . 
Along with it the reverse process of the demagnetization 



of almost homogeneous initial state on heating is possi- 
ble in the infinitesimal field. It is shown in Fig. [Dby the 
directed dashed lines. The presence of some small field 
here is necessary to remove the ferromagnetic instability 
of disordered states, without the field the inhomogeneous 
states will jump right to the ferro-state due to this insta- 
bility. 

The exhaustive study of both these processes was made 
in the graphite-intercalated magnet Cuq. 93(700.07(^^2 
having the phase sequence of path B^^ . The similarity of 
the experimental data to rh shown in Figs. I13[ [5lis quite 
impressive as authors has revealed the existence of nu- 
merous metastable states using multiple cooling- hea ting 
cycles, see the example in Fig. [T^l Also in Ref.|23j 
"thermal remagnetization" similar to that in Fig. [T51 
was observed in disordered ferromagnets Nd3{Fe,Ti)2g 
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FIG. 13: (Color online) Temperature dependencies of metastable states' magnetizations, susceptibility and heat capacity in a 
field h — a for j3 — 0.1 corresponding to the path B. 
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and Nds^Fe, Re)29- The depolarization and polariza- 
tion on heating were also seen in PLZT ferroelectric 
ceramics^iiS^. 

At last we should mention the example of the relaxor 
the thermodynamics of which is definitely distinct from 
that of the present model. It is PMN-PT with lead ti- 
tanate (FT) content between 0.06 and 0.2 in which the 
temperature evolution of hysteresis loop shows the grad- 
ual transition from the inclined to the rectangular form^ 
instead of the sharp one as in Fig. [HI This means that 
between glass and ferro-phase in this compound there 
are the sequence of mixed phases where the most stable 
states are those with partially ordered dipole moments. 
Probably the model with the three- and four - mode in- 
teractions can describe such intermediate mixed phases. 

V. DISCUSSION AND CONCLUSIONS 

Thus the present thermodynamic theory can describe 
in a unified manner a wealth of nonergodic phenomena 
in disordered ferromagnets and ferroelectrics. It demon- 
strates what the effects of numerous metastable states 
can actually be and how they can manifest themselves 
in various regimes of real-time quasi-static experiments. 
Also it shows what a form the results of rigorous statis- 
tical mechanics should have to describe the phase tran- 
sitions in crystals with random interactions. 

The important conclusion can be made on the rela- 
tions between the predictions of eqiulibrium thermody- 
namics valid in the infinite-time limit and finite-time ex- 
periments. Contrary to naive expectations that there are 
no such relations the present study shows that they can 



be revieled in finite-field quasi-static experiments. In par- 
ticular, they are straightfordly manifested in the change 
of hysteresis loop form closely related to the form of the 
static thermodynamic potential. Yet its form and ther- 
modynamic path for a given crystal can be obtained in 
some other sufficiently full set of finite-field quasi-static 
experiments. 'Sufficiently full' means here that a wide 
spectrum of inhomogeneous metastable states should be 
explored to determine the regions of their existence and 
their thermodynamic properties. The exellent examples 
of such studies providing the valuable data for the future 
theory are given in Refs. (isl. [2^. 

Apparantly the present phenomenology can be im- 
proved and expanded in many ways. Thus one may 
consider a multicomponent order parameter as, say, in 
Heisenberg ferromagnet or cubic ferroelectric and less 
sparse modes having three- and four-mode couplings in 
Landau's potential. Also one may include in it higher- 
order terms in m.i to expand the results to larger T and 
h regions. Eventually it may be worth to explore the 
models with the broken in some way permutation sym- 
metry. It seems also important to consider the role of 
fiuctuations (i. e. the modes which do not condense) 
in the present approach. The quantitative description of 
the nonergodic thermodynamics of PMN - type relaxors 
needs also the consideration of random-fields. 
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